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ABSTRACT 



A computer model is presented for determining the daily allo- 
cation of airstrikes for interdicting a lines-of-communication (LOC) 
network assuming an exponential damage function. The allocation of 
sorties to arcs in the network is based on the assumptions that (1) 
flow of supplies is restricted by network capacity and (2) parameter 
values for each arc are independent. Upper and lower flow bounds, 
cost per unit flow, maximum repair time and cost, and the vulner- 
ability parameter for each arc are required as input data. The 
model selects that arc to strike which maximizes the repair cost 
plus the product of the increase in minimum cost circulation flow 
and repair time. The procedure is programmed in daily cycles, 
allowing repair on interdicted arcs. A sample problem and all docu- 
mentation necessary for duplicating the computer program are given. 
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I . INTRODUCTION 



A. GENERAL 

The ability of a military force to conduct sustained combat 
operations any distance from its base of supply may be limited by 
the land lines-of-communication (LOC) network over which its 
supplies must flow. This is particularly true if air and sea 
resupply are not feasible, or the opposing force has superiority 
over these alternatives. The amount and type of material required 
to support a military force are determined by the nature of combat 
operations and the intensity of combat. A large conventional force 
would require more pounds-per-man-day than a small guerilla force 
capable of obtaining a portion of its requirements in the area in 
which it operates. Undoubtedly, however, a requirement will exist 
for supplies in excess of what can be obtained locally or what can 
be transported organically by the combat force. 

The land LOC network then becomes an important military target. 
Flow of military material over the network can be reduced by any of 
three alternatives: (1) attacking the transport capability, i.e., 

vehicles or rolling stock, (2) attacking manufacturing and supply 
centers, thereby reducing the amount of materiel available, or (3) 
attacking the road/rail LOC network over which the materiel must flow. 
Certainly, a combination of these alternatives would be most desirable. 
However, political considerations might preclude attacking manufactur- 
ing and supply centers; and terrain, weather, and lack of intell- 
igence information might limit the effectiveness of attempting 
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to attack convoys moving over the network. 

Such a scenario is particularly applicable to the present war 
in Southeast Asia and is a general characteristic of limited warfare. 
The North Vietnamese neither possess a sophisticated air resupply 
capability, nor can they openly move supplies by sea to support 
their forces operating in South Vietnam. They do, however, enjoy 
immunity from attack on all their facilities north of the 17th 
Parallel. Therefore, the supply of materiel to their forces oper- 
ating in the South is primarily limited by the transportation system 
over which this materiel must flow. 

The objective of an interdiction campaign by airs trike must 
then be to systematically reduce the capacity of road or rail 
segments which will ultimately restrict the amount of materiel that 
can flow through the network. 

B. AVAILABLE LITERATURE 

Several authors have published works on the problem of inter- 
dicting a lines-of-coramunication network. Durbin [3] has construct- 
ed a model for determining maximum cargo flow as a function of 
available transport vehicles. His method sequentially selects and 
interdicts the most vital segment in the network until flow is re- 
duced to zero or a specific number of segments have been destroyed. 
The model operates in daily cycles and allows repair on interdicted 
segments . 

Wollmer [11, 14, 15] has developed a procedure to determine 
the n most vital segments in a network and to evaluate the effect 
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of removing the most vital segment. In the latter work, road 
capacity as a function of the number of strikes is reduced (1) 
by a known quantity and then (2) by a quantity that is a random 
variable. 

Wollmer [12, 13] has also developed algorithms and a model 
for evaluating a targeting strategy whose aim is the greatest re- 
duction in network flow. The model is programmed in daily cycles 
with the user specifying the number of days and strikes. 

Mustin [9] uses the topological dual to represent a lines-of- 
communication network and assumes the decrease in capacity of a 
road or rail segment has a fixed deterministic rate which is linear 
between its upper and lower bounds. 

Nugent [10] also uses the topological dual for network rep- 
resentation, but assumes the reduction in capacity of a road or 
rail segment is exponential with decreasing marginal returns. His 
model is constrained by the number of available sorties and does 
not allow for the dynamics of the daily allocation problem. 

C. ALLOCATION OF EFFORT 

Many techniques have been devised to aid in solving the 
problem of allocating limited resources. Hancock [8] provides 
a reference on classical optimization techniques; Danskin [1, 2] 
considers the theory of maxima and minima and its application to 
weapons systems- allocation; and Koopman [6, 7] considers a simi- 
lar problem in the allocation of search and screening effort. 

The problem facing the strike planner is to determine what road 
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or rail segments in the network to attack and how many strikes to 
allocate to each selected segment. Numerous considerations must 
be weighed in order to make this determination: (1) the number of 
sorties available on a given day, (2) the probability that a strike 
against a particular segment is successful, (3) the expected cost 
to the enemy in time and effort to repair an interdicted segment, 
(4) the expected cost of aircraft lost in striking a defended seg- 
ment, and (5) the expected reduction in capacity if a strike is 
successful. 

The obvious procedure in making the allocation of W sorties 
is to attack those segments (at a level of n_^_. each) that maximize 
the ratio of benefits to cost. As will be shown later, Gibbs 1 
Lemma [2] provides the mathematical basis for determining this 
allocation. 



\ 
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II, OBJECTIVE AND SCOPE 



The objective of this paper is to present a model for determin- 
ing the daily allocation of airstrikes for interdicting a lines-of- 
communication network assuming an exponential damage function. The 
model seeks to allocate a maximum of N (specified) sorties per day 
to those arcs (segments) which yield the greatest reduction in flow 
capacity and the greatest time and cost to repair. Two options are 
provided that will allow the strike planner to (1) specify the 
number of strikes per day and the number of days in the interdiction 
campaign, or (2) select a percentage of the maximum possible flow 
to be reduced. Using the second option, the maximum available 
number of sorties per day and the number of campaign days are 
specified. 



Certain information is assumed to be available to the strike 
planner for each arc in the LOC network: 

(1) Upper flow bound - maximum capacity in units per day, 

(2) Lower flow bound - minimum capacity in units per day below 
which the flow cannot be reduced. 

(3) Maximum repair time - the time required to repair the 
maximum possible damage on an arc. 

(4) Maximum repair cost - the cost to repair the maximum 
possible damage on an arc. 

(5) Vulnerability parameter - the probability that a sortie is 
successful. The strike planner can manipulate this para- 
meter by providing his own definition of "successful." 



The vulnerability of any arc, , is defined as the upper flow 
bound minus the lower flow bound, i.e., U^-L„. A necessary assump- 
tion is that the vulnerability, repair cost, and repair time are in- 



dependent of the same parameters for any other arc. 
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By assuming an exponential damage function, the reduction in 



capacity d , for any arc (i,j) caused by n strikes is: 

d. . = V, .11 - exp (-VP . . n. .)] 
i j ij 13 13 

Arc repair time and cost are linear functions of the amount 
of damage caused by n. . strikes: 

13 

t.. = (d . . /V. . )T . . 

13 13 13 13 

k., = (d . . /V . )K. . 
ij 13 i 3 13 



( 1 ) 



( 2 ) 

(3) 



The amount of arc repair that can be completed in any one day 
is a function of the computed cumulative arc repair time and the 
total damage on the arc: 

r.. = (1/CRT..) D ±j (4) 

The value of allocating n _ strikes per day to arc (i,j) 



is determined by: 



Value. . = (F T - F) t. . + k. . (5) 

ij ij ij 

The model daily selects and interdicts those arcs that 
maximize (5) , performs repair as indicated by (4) , and steps to 
the next day. The procedure has been programmed in FORTRAN IV 
for use on the IBM 360/67 computer* 
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III. NETWORK REPRESENTATION OF LINES-OF-COMMUNICATION 



A, NETWORK STRUCTURE 

A transportation system, whether single or multimode, can be 
represented by a connected graph of nodes and directed arcs. For 
example, a highway system would consist of roads (arcs) and the 
intersection of roads (nodes). If a road is one-way, the beginning 
node of an arc is the point where traffic enters the segment and the 
ending node is the point where traffic leaves the segment. If a 
road is two-way, it is represented by two arcs, one in each direction. 

The LOC user attempts to send units of flow, expressible in 
vehicles per day, tons per day, etc., from the source to the sink. 

The source of flow could be a rear-area supply center and the sink 
could be the battlefront or a forward-area supply center. If the 
user is constrained by the capacity of the transportation system, 
it is assumed that he will attempt to maximize the total flow of 
cargo through the network at a minimum cost. 

Consider the example LOC networ' shown in Fig. 1. Suppose the 
user wishes to send flow from node AA to node ZZ. Node AA is the 
source and node ZZ is the sink. The route (AA,A,K,L,ZZ) is one of 
many possible routes over which the flow can be sent. Associated 
with each arc are certain parameters: lower bound on flow, upper 
bound on flow (flow capacity), and the cost per unit flow. For arc 
(K,L) the following parameters are given: 



Lower bound (L . . ) 0 

ij 

Upper bound (U^) 400 

Flow cost (C,^) 7 
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These parameters imply that between 0 and 400 units of flow per time 
unit can be sent over (K,L) at a cost of 7 per unit. This cost can 
be expressed in terras of dollars per ton of cargo, dollars per 
vehicle, manpower units, or any other definable cost. If the arc 
is used to its capacity, i.e., M^. = , it would cost 2800 per time 

unit to send flow over this arc. The total cost per time unit of 
sending 400 units of flow from node AA to node ZZ over route 
(AA,A,K,L,ZZ) would be: 

Z = 400(0) + 400(6) + 400(7) + 400(0) 

all (ij) 1J 1J 
in 

the route « 4200 

It is important to note that the flow out of a node equals the 
flow into a node. This conservation of flow implies there can be no 
storage at nodes. Therefore, for most general networks, if X units 
of flow are sent from the source to the sink, then any flow pattern 
that accomplishes this will have net flows of X units out of the 
source, X units into the sink, and zero on all other nodes. 

In actual LOC networks, there may be more than one source and/or 
sink. This can be handled, provided any source can supply any sink, 
by adding artificial nodes and arcs. This allows the user to con- 
vert the multi-source, multi-sink problem to a single-source, single- 
sink problem. The artificial arcs are used to connect the true 
sources (sinks) to an artificial super source (sink) . Note in 
Figure 1 that arcs (AA,A) , (AA,B), (AA,C), and (AA,D) can be con- 
sidered artificial arcs connecting the true sources A, B, C, and D 
to the single artificial super source AA. Similarly, arcs (L,ZZ), 
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(H,ZZ), (I,ZZ), and (J,ZZ) can be considered artificial arcs connect- 
ing the true sinks L, H, I, and J to the artificial super sink ZZ. 

The flow on an artificial arc, e.g., (AA,A), can be considered the 
rate that cargo leaves the true source A. The capacity for arti- 
ficial arcs is usually infinite, but could represent the maximum 
supply rate at a facility such as a supply center or warehouse. 

The lower flow bound and cost on artificial arcs are normally zero. 
Artificial arcs are not considered vulnerable to attack. 

Actual flow on an arc is designated M if the flow is from 

node i to node j and M.. if the flow is in the reverse direction. 

Ji 

Normally, flow is assumed to be from the source to the sink but 
can be in either direction on intermediate two-way arcs. 

, U _ , and represent the upper flow bound, the flow 
capacity at any time, and the lower flow bound, respectively, for 
any arc (i,j). The following relationship must hold for these 
parameters : 



0 < L. . < M. . < U. . < U. . 

- ij - ij - ij - ij 

The maximum possible flow through the network can be determined 
by using the maximal-flow minimal-cut theorem 14]. It states that 
the maximum flow in a network is equal to the minimum value of all 
cut-sets, where a cut-set is defined as a set of arcs that separate 
the source from the sink. The value of a cut-set is defined as the 
sum of the capacities of all arcs in the cut-set. In the example 
network, Figure 2, arcs (A,K) , (K,E) , (E,H) , (E,F), (C,F), (C,G), 
and (D,G) form the minimum cut-set with a value of 5000 units. 
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Figure 2. Minimum Cut-Set for Example Network 



R. CIRCULATION FLOWS IN NETWORKS 



As previously stated, if the LOC user is constrained by the 
capacity of the network, it is assumed he wishes to maximize the 
total flow through the network at a minimum cost. Such a problem 
is efficiently handled by the "Out-of-Kilter" algorithm developed 
by Fulkerson 14]. The network is restructured by adding a ,f universal M 
arc from the super sink to the super source as in Fig, 3, This 
artificial arc is assigned an upper flow bound greater than or equal 
to the maximum possible flow through the network (usually infinite), 
a lower bound of zero, and a unit flow cost less than the flow cost 
on any arc (usually minus infinity). Hence, the problem of maxi- 
mizing source-sink flow is converted to one of finding a minimum 
cost circulation flow by maximizing flow on the universal arc. 
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Figure 3, Example Network with Universal Arc 



IV. MODEL OPERATION 



A. MATHEMATICAL PRELIMINARIES 

Assume: (1) The set X* = (x* x* , x*) 

1 Z n 

n 

maximizes E f.(x.) 

i=l 1 1 

n 

subj . to E x. = X ;x. > 0 

i=l 1 1 ' 

(2) The are continuous and differentiable 

(3) f” (x.) < 0 

Note: These are necessary and sufficient conditions 

for a global maximum since the objective func- 
tion is concave and the constraints are convex. 

Then there exists a scalar s > 0 such that: 

(1) f.(x.) = s for x* > 0 

l i l 

(2) f](x.) < s for x* = 0 

ii i 



Note: The Lagrange Multiplier, s, can be interpreted 

as the implicit marginal price of allocating an 
additional unit of resource. 



The above is known as Gibbs’ Lemma and is discussed in more detail by 
Danskin [2]. It can be derived from the Kuhn-Tucker conditions or by 



elementary arguments based on the special structure of the problem. 



There are two important implications of this lemma: 



( 1 ) 

( 2 ) 



Marginal return must be above a threshold value, s, for a 
resource to be allocated 



t ^ 

For alternatives which receive allocation, f.(x.) = 

i.e., at the optimum allocation point, the 

ratio (marginal return/marginal cost) ^ is constant. 
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B. MODEL DESCRIPTION 



As stated earlier, there are numerous considerations to be 
weighed in allocating strikes to arcs in the LOC network. From a 
mathematical viewpoint, it is assumed that the effectiveness of the 
network to its user may be measured in terms of how to obtain maxi- 
mum flow at minimum cost. Therefore, the criterion to be used to 
allocate sorties is to minimize LOC effectiveness over time, taking 
into account the cost to repair damage on an arc. This is accomplish- 
ed by allocating strikes to the set of arcs that maximize the costs 
to the user. 

The effect of striking an arc is to reduce its capacity and 
increase the cost of a minimum cost circulation flow for a specified 
period of time. Therefore, in allocating N strikes, the selected 
arcs are the ones that maximize repair cost plus the product of the 
increase in minimum cost circulation flow and repair time. Specific- 
ally, if F is the cost of a minimum cost circulation flow before arc 
(i,j) is struck, F ? is the cost of a minimum cost circulation flow 

after n. . strikes on arc (i,j), and t . . and k. . are the repair time 
ij ij ij 

and repair cost, respectively, for damage caused by n_ strikes, 
then the arcs selected for attack will be the ones that maximize: 



Z (F* - F) t 
all (i,j) 




s.t. (1) Z n , . = N ; n. . > 0 
.n r-t *0 *3 ~ 



all (i,j) 



where : 
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F = Z M. .C. . - d..(C.. + C u ) ; C.. > 0 

all (i,j) 1J lj 1J 1J ^ 

C U « 0 



Note: C U is the unit flow cost on the universal 

arc* 



d ij = V 1 - V 3 

‘ll = <4 1J /V u > T « 

k. . = (d. ./V. .) K. . 
ij ij 13 13 

The objective function then becomes to maximize: 



.,Z,. .x-T..V. .(C. . + C U ) [1 - exp (-VP . . n. ,)] 2 
all (i,j) 1] i] i] v 13 13 



+ K. . [1 -exp(-VP. . n. . ) j 
13 13 13 

This can be shown to be a continuous concave function of n. . whose 

13 

second derivative is negative. The constraint is convex; therefore 
the above objective function becomes the f^(x^) of Gibbs' Lemma. 

To determine the allocation of W airstrikes, it is necessary to 
compute the value of the Lagrange Multiplier, s*. By Gibbs' Lemma, 

JL 

it is known that there exists an s >0 such that: 

n.. =0 for f!(X.) < s* 

ij 11 

n. , > 0 for f!(x ) > s* 

13 1 i 

An iterative method for this computation is developed by Nugent [10] 
The above procedure is used in Option (1) where the number of 
strikes is specified. For Option (2), a similar procedure is 
developed with the addition of another constraint. Recall that 
Option (2) requires a percentage by which the initial flow is to 
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be reduced. The additional constraint becomes: 

(2) Z M > (1 - P )m! 
all (i,j) 13 " 13 

I 

where: M. . is the initial network flow, 
ij 

p is the percentage by which the 
initial flow is to be reduced. 

Obviously, for the allocation of a maximum of N sorties, either 
constraint (1) or (2) will be active and the value of the Lagrange 
Multiplier, s*, can be computed using Nugent’s method. 

If the vulnerability parameter of arc (i,j) takes cognizance 
of the probability of losing an aircraft on any mission against 
that arc, then the model seeks to interdict that set of arcs where 
the ratio (marginal return/marginal cost)., is constant. 



C, COMPUTER SOLUTION PROCEDURE AND OUTPUT 



The computer solution procedure operates in daily cycles, the 
strike planner specifying by Option (1) or (2) the specified or 
maximum number of strikes per day, respectively, and the number of 
days in the interdiction campaign. The computer solution differs 
from the optimal model described in B., in that strikes are allocated 
one at a time. Each time, the arc selected to be struck is the one 
which maximizes : 



(F - F) t. . + k. . . 

ij iJ 

? 

The values F , F, t.., and k.. are computed for one strike instead 

ij 1J 

of for n„ strikes. This is a marginal allocation procedure, the 
solution having the property that if an arc is selected for inter- 
diction when W strikes are allocated, it will also be selected for 
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interdiction when more than W strikes are allocated. Since optimal 
allocation of strikes would not necessarily have this property, the 
selected arcs are optimal for one strike, but only approximate optim- 
ality for multiple strikes. 

After each strike, the following information is output: 

(1) Arc selected for strike 

(2) Arc capacity before and after strike 

(3) Actual arc flow before and after strike 

(4) Day arc repair is to be completed 

(5) Cumulative repair cost 

(6) Network flow after strike 
*(7) Network flow cost after strike 

(8) Number of strikes against that arc for that day 

*Note: Since the unit flow cost on the universal arc is an 

artificial cost used to convert a particular type of source- 

sink flow problem to one of circulation flow, it is not con- 
sidered in the printed output calculation. 

The end of a campaign day is determined by the option used. 

Using Option (1) , the day terminates when W sorties have been 
allocated. Using Option (2), the day terminates when either W 
strikes have been allocated or the total network flow has been re- 
duced to a specified percentage of the initial flow. The following 
information is output at the end of each day: 

(1) Arcs selected for strike 

(2) Number of strikes allocated against arc (i,j) for 
that day 

(3) Number of strikes allocated against arc (i,j) for 
the campaign 

(4) Total number of strikes allocated that day 

(5) Total number of strikes allocated for the campaign 

At the start of each day, all arc repair scheduled for that 
day has been completed. The resulting flow pattern is computed and 
arc flow, maximum network flow, and network flow cost are output. 
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v. 



SAMPLE PROBLEM 



The network shown in Fig. 4 is used as a sample interdiction 
problem. It is the same network used previously but with the 
addition of maximum repair time and maximum repair cost parameters 
for each arc. Figure 5 depicts the same network, but node positions 
in relation to each other have been changed to indicate a more 
realistic lines-of-communi cation network. Artificial arcs from 
the super source are used to indicate that flow in the enemy sanctu- 
ary area is invulnerable to attack. Nodes A, B, C, and D are points 
on supply routes where arc capacity is subject to interdiction. The 
enemy objective is to move supplies from his sanctuary area to for- 
ward supply centers located at nodes L, H, I, and J. Nodes K, E, F, 
and G are route intersections. Figure 5 also gives the actual 
direction and magnitude of flow before the interdiction campaign 
begins . 

The objective of the strike planner is to reduce the initial 
maximum network flow by 60 percent and thereafter, maintain this 
reduced flow for a campaign of 5 days duration. He has a maximum 
of 50 aircraft available each day to fly interdiction missions. 

The computer solution to this problem is given in Appendix D. 
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Figure 4. Network for Sample Problem 
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VI. SUMMARY AND EXTENSIONS 



A computer model has been presented for determining the daily 
allocation of airstrikes for interdicting a land lines-of- 
coramunication (LOC) network, assuming an exponential damage function. 
The allocation of sorties to arcs in the network is based on the 
assumptions that (1) flow of supplies is limited by network capacity, 
and (2) parameter values for each arc are independent. The model 
selects that arc to strike which maximizes the repair cost plus the 
product of the increase in minimum cost circulation flow and repair 
time. The procedure operates in daily cycles, performing repair on 
arcs that have been struck. 

Two options are provided that allow the strike planner to 
(1) specify the number of strikes per day and the number of days 
in the campaign, or (2) prescribe a percentage by which the initial 
flow is to be reduced, the maximum number of strikes per day, and 
the number of campaign days. The upper and lower flow bounds, 
cost per unit flow, maximum repair time and cost, and the vulner- 
ability parameter for each arc in the network are required as in- 
put data. 

The model has been programmed in FORTRAN IV. for use on the 
IBM 360/67 computer. Detailed information for each day in the 
campaign is listed in the computer output to include arcs selected 
for interdiction, network flow and network flow cost at the end 
of each day, and the number of sorties allocated to selected arcs. 
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There are almost an endless number of extensions that can be 



considered in the allocation of airstrikes to a LOC network. In 

the present model, it is assumed that only one type of aircraft is 

available to attack the network, i.e., one weapons system against 

y targets, A more practical and interesting approach is to assume 

that more than one type of aircraft is available for interdiction, 

each with a different kill probability against different types of 

targets, i.e,, X weapons systems against y targets. The damage on 

an arc caused by n. . strikes then becomes: 

ij 



d. . 
iJ 



il 



[1 - it exp (-VP. 






Another modification would be to vary the damage function it- 
self. The exponential was used in this model because it exhibits 
decreasing marginal returns to scale and simplifies computational 
procedure; however, other such functions could be utilized. The 
exponential is a continuous approximation, for n _ large and VP_ 
small, of the function: 



V 1 - (1 - y V" ij;i 

A contrast could be made of the damage caused by both distributions 
with actual damage data. Holliday 1 5] has done* research in the 
field of damage assessment for the Advanced Research Projects 
Agency and his work might be used for such a contrast. 
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APPENDIX A 



COMPUTER PROGRAM INPUT DATA 



Card 
no . 


Column 
no . 


Format 


Description 




1 


1-80 


A 


Program title, (PRONAM) 




2 


1-5 


15 


Maximum number of strikes 
allocated each day, (N) 


that can be 




6-10 


15 


The number of days in the 
(NUMDAY) 


campaign. 




11-15 


15 


The number of invulnerable or artificial 
arcs in the network. Include all arcs 



considered invulnerable to airstrike, 
artificial arcs connecting the dummy 
source to the true source(s), and arti- 
ficial arcs connecting the true sink(s) 
to the dummy sink. (INVL) 





16-20 


15 


If equal to 1, specifies one-way flow on 
some arcs in the network. If equal to 
2, specifies two-Way flow on all arcs in 
the network. When arcs are two-way, a 
strike against one implies a strike 
against the other. (X2WAY) 




22-26 


F4.2 


Specifies the percentage of flow to be 
interdicted. If greater than 1.0, a 
fixed number of strikes per day (N) will 
be targeted. (REDU) 


Arc 


1-4 


A4 


Name of arc. If two-way flow, must be 


Data 






the same for both arcs. (ARNAM) 


Cards 


6-7 


A2 


Beginning node. (FROM) 




9-10 


A2 


Ending node. (TO) * 




12-20 


19 


Lower bound on arc flow, (LBND) 




22-30 


19 


Arc capacity (upper bound on arc flow) • 
(MCTY) 




32-40 


19 


Cost of one unit of flow on the arc. 



(JFCST) 



30 



42-50 



19 



52-60 19 



66-70 F5.3 

Last 

Card 1-2 A4 



Maximum arc repair cost, i.e., the cost 
of repairing the arc from its lower bound 
capacity to its maximum capacity. (KMAXCT) 

Maximum arc repair time, i.e., the time 
required to repair the arc from its lower 
bound capacity to its maximum capacity. 
(MRTM) 

Arc vulnerability parameter. (VP) 



The word ’GO' . 



Note: Arc data cards must be placed in the following order: 

1. Universal arc. 

2. Invulnerable and artifical arcs. 

3. Vulnerable arcs. 
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APPENDIX B 





GLOSSARY OF PROGRAM TERMS 


ARNAM(k) 


The designation of arc ij . (input data) 


COSTMN(k) 


Real number equivalent of MNCOST. 


CRC 


Cumulative repair cost for all interdicted arcs. 


CRT (k) 


Cumulative repair time for an interdicted arc. 


CSTMNT 


Real number equivalent of MNTCST. 


DAMAGE (k) 


Total reduction in capacity of arc ij caused by one 
or more strikes against it. 


DAY 


Counter for number of days up to NUMDAY. 


FLOW(k) 


Amount of flow on arc i j . 


FROM(k) 


Location of beginning node of arc ij . (input data) 


GO 


Flag used to denote last data card has been read. 


IKEEP1 


Used to save capacity of arc ij . 


IKPOP2 


Used to save capacity of the optimal arc to strike. 


IKPOP1 


If OPTARC is two-way, saves the capacity of the 
reverse direction arc. 


IMNCST 


Cost of the initial flow in the network. 


IMNCUT 


Cost of flow in the network after the optimal arc 
has been struck. 


INFEAS 


A flag used in the Out-of-Kilter subroutine to 
denote infeasible flow. 


INVL 


The number of invulnerable or artificial arcs in 
the network, (input data) 


IOPARC 


If an arc is two-way, specifies the reverse 
direction of the best arc to strike. 


IRQFLO 


Integer equivalent of REQFLO. 
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IS(k) 


Counter for arcs that have been struck one or more 
times , 


ISWTCH 


Flag used to denote a two-way arc. 


IT(k) 


Counter for total strikes on arc ij during the 
campaign. 


IV (k) 


Counter for total strikes on arc ij during one day 
of the campaign. 


JFCST(k) 


Cost of one unit of flow on arc ij . (input data) 


JFLOW(k) 


Saves flow on OPTARC before strike, (for output) 


J2WAY 


Same as ISWTCH. 


KFLOI 


Saves flow on OPTARC after strike, (for output) 


LBND (k) 


Lower flow bound on arc ij , (input data) 


LINE 


Line counter used in the output subroutine. 


MCAP (k) 


Capacity of arc ij before next strike. 


MCAPU 


Saves the capacity of OPTARC before strike. 


MCTY (k) 


Flow capacity of arc i j . (input data) 


MNCOST 


Network flow cost before next strike. 


MNTCST 


Network flow cost after arc ij has been struck. 


MRTH(k) 


Maximum repair time for arc ij . (input data) 


N 


Maximum number of strikes that can be allocated 
during any one day. (input data) 


NEWCAP (k) 


Capacity of arc ij after strike. 


NN(k) 


Unique node number assigned by Subroutine Number. 


NR 


Counter for total strikes on all arcs during any 
one day. 


NRC 


Counter for total campaign strikes. 


NUMARC 


Maximum number of arcs in the network. 


NUMDAY 


Number of days in the campaign. (input data) 
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NUMNDE 

OPTARC 

REDU 

REPDAY(k) 

RMAXCT(k) 

RPMCST (k) 

TME(k) 

TO(k) 

TOTCST 

TOTDST 

VALUE (k) 

VP(k) 

VULCAP (k) 

XINCRE(k) 

X2WAY 



Maximum number of nodes in the network. 

The optimal arc to strike. 

The percentage by which initial flow is to be 
reduced. (input data) 

Day on which arc ij is to be repaired. 

Maximum repair cost for arc ij . (input data) 

Marginal repair cost for arc i j . 

Marginal repair time for arc ij . 

Location of ending node of arc ij . (input data) 

Counter for cumulative campaign strikes. 

Counter for cumulative day strikes. 

The value of striking arc ij , 

Vulnerability parameter of arc ij . (input data) 

The vulnerable capacity of arc ij . 

Marginal reduction in capacity on arc ij caused by 
an additional strike. 

Input parameter that specifies whether or not all 
arcs are two-way. 
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A P P H N D I X C 



PROGRAM LISTING 



DIMENSION XINCRFU50) 

COMMON /CAT 01 / I NVL , DA Y , NUMDAY , NUMARC , OPT ARC , L I NE , INFcAS 
COMMON/ DAT02/S»NUMNDE* M CAPU*KFLCI , NR , N 
COMMON /DA TC3/NRC , I W NCST , I MNCUT »X2WAY » M NCOST f MNTCST 
COMMON /CAT 04/REPDA Y ( 13 0 ) ,ARNAM( 150 ) , FROM I 1 50 ) , TO I L5 0 ) 
C0MM0N/DAT05/RMAXCT( 150) , R P MC S T ( 15 0 ) , P I ( I 5 0 ) 

COMMON /CAT 06 /LB NCI 150 ) ♦ NE WCP ( 1 50 ) , MC AP ( 150 ) , I S ( 1 50 ) 
COMMON/DAT 07 /MRTMI 15 0 ) , NN( 1 5 0 ) , MCT Y ( 1 50 ) , J F CS T ( 1 50 > 
C0MM0N/CAT08/FL0WI 150) ,NL( 150) 

COM MON/DAT 09 /VALUEI150) , PRGNAMI 25 ) , JFL OW ( 15 0 ) 

COMMON /DAT 10 /REQFLO ,REDIJ, J2WY, I SWTCH ,CPC 

COMMON /DAT 11/ I OP ARC, TOTDSTI 150) , TOTCST ( 150 ) , I T ( 1 5 0 ) 

COMMON /CAT 12 /I VI 150) , DAMAGE (15 0) , VP I 15 0) 

C0MM0N/DAT13/ I RQFL C , VULC AP ( 150) , CRT ( 150) 

COMMON /D AT 1A/TMEI 150) , VULC TY( 150) 

INTEGER TOTDST, TOTCST, CAY, ARNAM,X2 WAY, FROM, TO 
INTEGER FLOW, PI , OP TARC , REP DAY , S , RMAXCT ,CRC 
DATA NAME/ 1 GO ' / 



C 



C 



C 



C 

C 



OPTFLO = 0. 

REQFLO = 0. 

IRQFLO = 0 
1 = 0 

READ PROGRAM TITLE AND PROGRAM PARAMETERS 

READ(5,5C5)(PRONAM(KK) ,KK=1,20) 

READ I 5 , 500 N, NUMDAY, I NVL ,X2WAY , REDU 
WRITE! 6,5CA) 

S = INVL + 2 

READ ARC PARAMETERS 

10 1=1+1 

READ(5,51C)ARNAM( I ),FROM(I ) , TO ( I ) , LBND ( I ) , MCTY ( I ) 

1 , JFCST I I ) , RMAXC T ( I ) , MRTMI I ) , VP I I ) 

IFIARNAM(I) .EQ. NAME) GO TO 20 

DETERMINE PARAMETERS FCR FIRST STRIKE 

VULC AP I I ) = MCTYII) - LBND ( I ) 

MC AP I I ) = MCTY I I ) 

CAMAGE I I ) = VULCAPII )*(1.0-EXP(-VP(I)) ) 

NEWCPII) = FLOAT (MCTY! I) ) -DAMAGE ( I ) + .5 
RPMCSTII) = (FLOAT! R^AXCTI I ) J/VULCAPI I ) )*DAMAGE( I ) 
TMEII) = I FLOAT! MRTMI I ) ) / V'JLCAPI I) )*CAMAGE( I ) 

GO TO 10 
2 C CONTINLE 

NUMARC = 1-1 

CALL NUMBER SUBROUTINE TO DESCRIBE NETWORK 

CALL NUMBER 

INITIALIZE VARIABLES 

DO 30 K = 1, NUMARC 
ISIK) = 0 
ITIK ) = 0 
I V I K ) = 0 
FLOWIK) = 0 
NLIK) = 0 
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c 



c 



c 



c 



c 



c 



c 



c 

c 

c 



REPDAY(K) = 0 
TOTCST(K) = 0 
TOTDST(K) = 0 
CRT ( K ) = 0. 

VALUE ( K ) =C . 

VULCTY(K) = 0. 

XINCRE(K) = 0. 

30 CONTINUE 

SOLVE UNI NTERDICTfcD NETWORK 

CALL OCK( NUMNDS, NUMARC, FROM, TO, JFCST , MCAP , L RND , FLOW , 
1PI ,NL, INFEAS) 

I F { INFEAS. EO.l ) GO TO 70 
OPTFLO = FLOAT ( F LOW ( 1 ) ) 

RtQFLO = GPTFLO+l l.-REDU) 

IRQFLG = REOFLO 

PRINT PROGRAM TITLE AND MODE OF OPERATION 

CALL OUTPUT! 1 J 
CAY = C 
CALL OUTPUT ( 6 ) 

FIND MAX FLOW, FLOW COST, AND PRINT RESULTS 

CALL FLOCSTIFLOW, JFCST , NUMARC , IMNCST ,INVL, 1) 

CALL OUT PUT ( 7 ) 

DAY = 1 
NR = 0 
NRC = 0 
CRC = 0 

INTERDICTION OF NETWORK BEGINS HERE 

40 NR = NR + 1 
NRC = NRC + 1 
INFEAS = 0 
DO 45 I = S , NUMARC 
JFLOWl I ) = FLOW! I ) 

45 CONTINUE 

CALL STKARC TO DETERMINE BEST ARC TO STRIKE 
CALL STKARC 

IF ( INFEAS. EO. 1 ) GO TO 70 

OPTARC IS THE BEST ARC 

MCAP U = MCAP( OPT ARC ) 

KFLOI = F LCW ( OPT ARC ) 

ISAVE = OPTARC 
K = IOPARC 

STORE COST OF REPAIR AND COST OF STRIKE 

CRC = CRC + RPMCST (OPTARC) + .5 
50 CONTINUE 

CRT ( K I = CRT(K) + TMF(K) 

REPDAY(K) = DAY + CRT ( K ) ♦ .5 

COUNTER FOR ARCS THAT HAVE BEEN STRUCK 

IS ( K ) = I S ( K ) + 1 

COUNTER FOR DAY STRIKES 

I V ( K ) = I V ( K ) + 1 

COUNTER FOR CAMPAIGN STRIKES 

I T ( K ) = I T ( K ) + 1 
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ooooo o no 



TOTDST(K) = I V ( K ) 

TOTCST(K) = IKK) 

DETERMINE PARAMETERS FOR NEXT STRIKE 

J = I S ( K ) + 1 
MCAP(K) = NEWCP(K) 

VULCTY(K)=MCAP( K) - LBND(K) 
XINCRE(K)=VULCTY(K)+(1.-EXP(-VP(K) ) ) 

NEWCP( K)=FLCAT (MCAF(K) )-XINCRE(K) + .5 
DAMAGE ( K ) = MCTY(K) - NEWCP(K) 

RPMCST(K) = (FLOAT (RMAXCT(K) ) / VULC A P (K ) ) *X I NCR E ( K ) 
TME ( K ) =( FLOAT ( MRTM(K) ) / VULCAP ( K ) )* X INCR£(K ) 

IF ARC IS 2-WAY, GO BACK AND COMPUTE 
PARAMETERS FCR RETURN ARC 

IF ( ISWTCH.EO. 1) GO TO 60 
K = ISAVE 
ISWTCH = 1 
GO TO 5C 
60 CONTINUE 

PRINT STRIKE RESULTS 

CALL OUTPUT ( 2 ) 

7 C CONTINUE 

IFUNFEAS .EO. 1) GO TO 130 
I F( FLOW ( 1 ) .LT. LBND(l)) GO TO 130 
IF ( FLOW ( 1 ) .EQ. 0) GO TO 140 

IF REDL IS < 1, NUMBER OF STRIKES IS NOT 
CONSTANT. IF NUMBER OF STRIKES = NUMBER 
SPECIFIED IN INPUT OR IF CAPACITY 
HAS BEEN REDUCED TO SPECIFIED AMOUNT, 

DAY IS DONE. 

IF(REDU.LE.1.0)G0 TO 80 
I F ( NR .EO. N) GO TO 90 
GO TC 40 
80 CONTINUE 

IF(FLOW( D.LE.REOFLOJGO TO 90 
I F ( NR . EQ . N )G0 TO 90 
GO TO 40 

90 CALL OUTPUT ( 5 ) 

IF ( DAY .EO. NUMDAY) GO TO 150 
LINE = 48 

RE-INITIALIZE PARAMETERS FOR NEXT DAY 

100 DAY = DAY + 1 
NR =0 

DO 110 M = S ,NUMARC 
I V ( M ) = 0 
TOTDST(M) = 0 
110 CONTINUE 

PERFORM REPAIR CN ARCS 

DO 120 M = S ,NUMARC 
I F ( I S ( M ) .EQ. 0) GC TO 120 

MCAP ( M ) = FLCAT ( MCAP(M) ) + (l. /CRT ( M ) ) *DAM AGE ( M ) 
VULCTYt M)=MCAP( M) - LBND(M) 
XINCRE(M)=VULCTY(M)*(1.-EXP(-VP(M) ) ) 

NEWCP< M)=FLOAT ( MCAP< M ) ) - XINCPE(M) + .5 
DAMAGE(M)=MCTY(M) - NEWCP(M) 

RPMCST(M) = (FLOAT (RMAXCT(M) )/VULCAP(M) )*XINCRE(M) 
TME ( M)=( FLOAT ( MRTM ( M ) ) /VUL C AP ( M ) )*X I NCRE ( M ) 

CRT(M) = CRT ( M ) - 1. 

120 CONTINUE 
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C SOLVE NETWORK FOR NEXT DAY 

CALL OOK( NUMNDE,NURARCt FROM, TO, JFCST,MCAP, LB NO, FLOW, 
1 P I » NL » INFEAS) 

I F { INFEAS .EO. 1 ) GO TC 70 

C FIND FLCW COST 

CALL FLCCST( FLOW , JFCST, NUMARCt IMNCST, I NVL » 1 ) 

C PRINT NETWORK FLCW PATTERN AND ITS COST 

CALL OLTPLiT ( 6) 

CALL 0 U T P U T ( 7 ) 

IF(FL0W(1J .GT.REQFLOIGC TO 40 
GO TO ICO 
1 3 C CALL OUTPUT ( 3 ) 

GO TO 90 

140 CALL OUTPUT ( 4) 

GO TO 90 
150 CONTINUE 
CALL EXIT 



500 FORMAT (4I5flX,F4«2) 

5C4 FORMAT ( 1H1,1.0X) 

505 FORMAT ( 20A4 ) 

510 FORMAT (A4 t 2(lX t A2) , 5 ( 1 X , 1 9 ) , 5X , F5. 3 ) 
END 
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on no o on noon non noon no 



00K IS THE FULKERSON OUT OF KILTER ALGORITHM 



SUBROUTINE 00K( NODES , ARCS , I , J ,CCST , H I , LO , F LOW , P l 
1 , NA , INFEAS > 

DIMENSION I( 150) , J ( 1 50 ) , COST ( 150 ) , HI (150) , 

1L0 ( 150 ) ,FL0W( 150) 

DIMENSION PI ( 150) ,NA( 150) 

INTEGER NODES, ARCS, I , J, C OS T, HI, LC, FLOW, PI, NA, INFEAS 
INTEGER A,AA,N,SRC , SNK , DEL , INF,C»ACK , CCK , N1 ,N?, 
1INC, LABEL 

DEFINITION OF CALLING SEQUENCE 

NAME USE 

NODES NUMBER OF NODES 

ARCS NUMBER OF ARCS 

I LIST OF FROM (BEGINNING) NODES 

J LIST OF TO (ENDING) NOCES 

COST UNIT COST OF FLOW CN ARCS 

HI UPPER BOUNDS FOR ARCS 

LO LOWER BOUNDS FOR ARCS 

FLOW AMOUNT OF FLOW IN ARCS 

PI NODE PRICES 

NA SCRATCH LIST FOR NODES LABELING 

INFEAS FLAG DENOTING THE CONDITION OF OUTPUT 



INF = 
AOK = 



2147483647 

C 



LOOK FCR AN OUT OF KILTER ARC 
AA=1 

100 N 1 = I ( AA) 

N2 = J ( A A ) 

C = COST ( AA ) + PI(Nl) - P I ( N2 ) 

40 IF(FLOMAA) .LT.LO( AA) .OR. ( C. LT . 0 . AND .FLOW ( AA ) . LT . 
1HKAA) ) ) GO TO 50 

IF(FLOMAA).GT.HI(AA).CR.(C.GT.O.AND.FLOW(AA).GT. 
1L0(AA) ) ) GO TO 60 
90 AA = AA + 1 

IF(AA.LE.ARCS)GO TO 100 

NO GUT OF KILTER ARCS LEFT 

INFEAS = 0 
RETURN 

OUT OF KILTER ARC FOUND 



50 


SRC 




J ( AA ) 




SNK 




I ( AA) 




LABEL 1 


= +AA 




GO 


TO 


2C0 


60 


SRC 


= 


I ( AA ) 




SNK 


= 


J ( AA ) 




LABEL 


= -AA 



SAVE LABELS IF LAST OPERATION WAS INCREASING NODE 
PRICES CN THIS ARC 



200 



IF( AA.EQ.ACK.AND.NA(SRC) .NE.O) 
DO 201 N =1, NODES 
NA ( N ) = C 
201 CONTINUE 
AOK = AA 
205 COK = C 

NA(SRC) = LABEL 



GO TO 205 
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on 



C LABEL 

2 1 C LABEL = 0 

00 250 A = 1 » ARC S 
N1 = 1(A) 

IF(Nl.LT.O) GO TO 250 
N2 a J ( A ) 

IF(NA(Nl).EQ.O.AND»NA(N2).EQ.O) GO TO 250 
IF(NA(N1).NE.0.AND.NA(N2) .NE.OJGO TO 245 
C = COST ( A ) + PI INI ) -PI (N2) 

IF ( NA( M I .EQ.O) GO TO 220 

I F ( FLOW ( A ) .GE.HI (A ) . CR . ( FLOW< A).GE.LO( A) .AND.C.GT .0) ) 
1G0 TO 245 
NA(N2) = A 
GO TO 240 

220 IFiFLOM A) . LE . LO ( A ) . OR . ( FLOW! AI.LE.HK A) .AND.C.LT.O) ) 
1G0 TO 245 
NA(N1) = -A 
240 LABEL = 1 

C NOOE LABELLED, TEST FCR BREAKTHRU 

IF(NA( SNK) .NE.O) GO TO 260 
245 1(A) = -N1 
250 CONTINUE 

GO BACK AND DO MORE LABELLING IF SOME NODE WAS 
LABELLED CN LAST PASS 

IF(LABEL.NE.O) GO TO 210 

C RESTORE POSITIVE SIGNS TO FIRST NODE LIST 

260 DO 270 A = 1 , ARC S 
1(A) = IABS(I(A)) 

27C CONTINUE 

C IF NO LABELLING DONE ON LAST PASS, GO TO 

C INCREASE PIE 

I F ( LABEL. EC. 0) GO TO 400 

C BREAKTHRU, FIND THE INCREMENT 

3CC INC = INF 

C FOLLOW PATH BACK FROM SOURCE 

N * SRC 

310 A = I A B S ( NA( N) ) 

IFINA(N).LT.O) GO TO 315 
N 2 = I (A) 

C = COST(A) - PI (N) + PI (N2) 

I F ( C.GT.O ) INC = HI N0( I NC ,L0( A) - FLOW(A)) 
IF(C.LE.O)INC = MINO(INC,HI (A) -FLCW(A)) 

GO TO 340 
315 N2 = J ( A ) 

C = COST( A) + PI (N) -PI (N2) 

320 IF(C.LT.O) INC = M I NO ( I NC , F LOW ( A ) — H 1(A)) 

IF(C.GE.O)INC = MINO( INC,FLOW( A) -LO ( A ) ) 

340 N = N2 

IF ( N .NE .SRC )G0 TO 310 

C INCREMENTS 

350 A = I A e S ( NA ( N ) ) 

I F ( NA < N ) .LT.O)GO TO 36C 
FLOW ( A ) = FLOW ( A ) + INC 
N = 1(A) 

GO TO 370 

3 6 C FLOW ( A ) = FLOW(A) -INC 
N = J ( A ) 
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3 7 C IF ( N .N E .SRC ) GO TO 350 

C FLOW INCREMENTED, RETURN TO KILTER TEST 

NA(N) = 0 
GO TO IOC 

C CHANGE PI 

ACC DEL = INF 

C FIND INCREMENT 

DO 420 A = 1 , A RC S 
N 1 = 1(A) 

N2 = J ( A ) 

IF(NA(N1)«EQ.0.AND.NA(N2).EQ.0)G0 TO 420 
I F ( NA ( N1 ) .NE.O.AND.NA( N2 ) . NE . 0 ) GO TO 420 
C = COST ( A ) + PI (Nl) - PI(N2) 

IF(NA(N2) .EQ.O.AND.FLQW(A).LT.HKA) ) DEL =MINO(DEL,C) 
IF(NA( N2) .NE.O.AND.FLOW( A) .GT.LO(A) ) DEL =M I NO ( DEL , -C I 
420 CONTINUE 

IF(DEL .NE.INF) GO TO 430 

IF(FLOW(AA).EQ.LO( AA ).CR.FLOW(AA). EQ . H I ( AA ) ) GO TO 425 

C INFEASIBLE SOLUTION 

INFEAS = 1 
RETURN 

C INCREASE PI 

425 DEL = lABS(COK) 

430 DO 450 N =1, NODES 

IF ( NA ( M • EC* 0 ) P I ( N ) = PI(N) + DEL 
450 CONTINUE 

C GO BACK TO KILTER TEST 

GO TO 100 
END 
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SUBROUTINE FLOCSTI FL CW , COS T , NUMARC ,MINCT, INVL , JUNIV) 
DIMENSION FLOW! 150) ,C0ST(150) 

INTEGER FL CW » COST 

MINCT = 0 
IUNV = INVL 

C IF JUNIV = 0, USE CCST OF UNIVERSAL ARC 

IF ( J UN I V .EO. 0) IUNV = I 
DO 10 I = I UNV » NUM ARC 
10 MINCT = MINCT + FLOW ( I ) *COST I I ) 

RETURN 

END 



SUBROUTINE NUMBER 

COMMON /D ATOl /INVL, DA Y.NUMD AY, NUMARC, OPT ARC, LINE, INFEAS 
COMMON/DAT 02/S, NUMNQE , MCAP U , KFL C I , NR , N 
C0MM0N/CATC4/REPDAYI ISO) .ARNAMI150) , FROM! 1 50 ) , TO I 150 ) 
COMMON /DAT07/MRTMI 1 50) , NNI 1 50 ) , MCT Y ( 150 ) , J FCST I 150 ) 
INTEGER TOTDST,TOTCST,DAY, ARNAM,X2WAY, FROM, TO 

C NUMBER THE NODES 

NUMNDE = 0 

DO 30 L = l, NUMARC 

NN ( NUMNCE + 1) =FRCM(L) 

FROM ( L ) = NODNUMIFRONIL) ,NN) 

NUMNDE = MAXO(NON,NDt»FRCM(L) ) 

30 CONTINUE 

DO 40 L *1, NUMARC 
NN( NUMNDE + 1 ) = TC( L) 

TOIL) = NODNUM (TOIL) ,NN) 

NUMNDE = MAXOI NUMNDE , TO! L) ) 

40 CONTINUE 

C CHECK FOR DEAD NODES 

DO 48 I =1, NUMNDE 
DO 42 L =1, NUMARC 
IF I FROM I L ) .EO. 1 ) GO TO 44 
42 CONTINUE 

WRITE! 6 , 522 ) NNI I ) 

44 CONTINUE 

DO 46 L =1, NUMARC 
IFITOILJ.EC.DGO TO 48 
46 CONTINUE 

WRITE! 6,530)NN( I ) 

48 CONTINUE 
RETURN 

522 FORMAT 1 1H1 , 22HN0 ARC BEGINS AT NODE , A4/1H1) 

530 FORMAT I 1H 1 , 22HN0 ARC ENDS AT NODE , A4/1H1) 

END 



FUNCTION NODNUMI 1 1 I ,NN) 

DIMENSION NN I 1 ) 

COMMON /CAT 02/S, NUM NDE,MCAPU,KFLCI , NR ,N 

NODNUM = NUMNDE + 1 
IF I NUMNDE . 60. 0 ) RETURN 
CO 10 I=i, NUMNDE 
IF! I 1 1 .EQ.NNI I ) ) NCDNUM = I 
1C CONTINUE 
RETURN 
END 



42 



on no 



SUBROUTINE STKARC 

C0MMGN/DAT01/INVL , DAY , NUMDAY , NUMARC , GPT ARC , L I NE .INFEAS 
COMMON/ CATO 2 /S , NUMNDE , MC A PU , KF L C I , NR , N 
COMMON /CAT03/NRC, I MNCST , I MNCUT , X2WAY , MNCOS T , MN TCS T 
COMMON / C A T 04/R fcP DA Y ( 150) ,ARNAM( 150 ) , FROM ( 1 5 0 > , TO ( 150 ) 
COMMON /DA T 05/RMAXCT ( 150) ,RPMCST(150).PI(150) 

COM MON /DA T 06 /L BND( 150) ,NEWCP( 150) ,MCAP( 150 > , I S ( 130) 
COMMON /DAT 07 /MR TM( 150),NN(150), MCTY ( 150) , JFCST ( 150) 
COMMON /DA TCS/FLOM 15C),NL(150) 

CO MM ON/CAT 09/VALUE( 150) ,PRONAM(25) ,JFL0W(150) 

COMMON /DAT 10 /REO FLO ,REDU, J2WY, I SWTCH.CPC 

COM MON/DAT 11/ IOPARCtTOTDST ( 150 > , TOTCST ( 1 50 ) , I T ( 150 ) 

COMMON /DAT 12/1 V( 15 0) , DAMAGE ( 150) , VP ( 150) 

COMMON /DAT 13/ I ROFL G » VULC AP ( 150) ♦ CRT ( 150) 
C0MMGN/DAT14/TME(150) , VULCTY( 150) 

INTEGER TOTDST , TOTCST , DAY, ARNAM,X2 WAY, FROM, TO 
INTEGER FLOW, PI , CP T ARC , REP DA Y, S , RMAXCT , CRC 

C INITIALIZE PARAMETERS 

MNTCST = 0 
MNCOST = C 
CSTMNT=0. 

COSTMN=C . 

UPVAL = -10000000. 

C SOLVE FCR LAST FLOW PATTERN 

CALL OCK( NUMNDE, NUMARC,FRGM, TO, JFCST, MCAP, LBND, FLOW, 
1PI ,NL, INFEAS) 

IF ( INFEAS. EQ. 1 JRETURN 

CALL FLCCSTI FLOW, JFCST, NUMARC, MNCOST, I NVL,0) 

DO 100 I = S , NUMARC 
I F ( F L 0 W ( I ) .EO. 0 ) GO TO 100 
IKEEP1=MCAP( I ) 

MCAP ( I ) =NE WC P ( I ) 

SOLVE INTERDICTED FLCW PATTERN FOR 
SELECTED ARCS AND FIND FLOW COST 

CALL OCK( NUMNDE , NUMARC, FROM, TO, JFCST , MCAP, LBND, FLOW, 

IP I , NL , INFEAS) 

I F ( INFEAS. EQ. 1JRETURN 

CALL FLCCSTI FL OW , JFCST , NUMARC , MNTC ST , I NVL , 0 ) 

FIND VALUE OF STRIKING ARCS AND 
SAVE MAXIMUM VALUE 

MCAP ( I ) = IKEEP1 
CSTMNT= MNTCST 
COSTMN=MNCOST 

VALUE ( 1 ) = (CSTMNT-COSTMN) *TME( I ) +RP MCST ( I ) 

IF(VALLt( I).LT .UPVALJGC TO 90 
LP VAL = VALUE ( I ) 

C OPTARC IS THE BEST ARC TO STRIKE 

CPTARC= I 
90 CONTINUE 
100 CONTINUE 

DO 52 I = S, NUMARC 

IF (FROM (OPTARC) . EO . T0( I ) . A ND . TO ( OP T ARC ) . EO . FROM ( I ) 

1 .ANO.ARNAM(OPTARC) .EO. ARNAM( I ) ) GO TO 54 

52 CONTINUE 

53 CONTINUE 
J2WY = 1 

C IOP ARC IS THE REVERSE DIRECTION ARC 

IOPARC = CPTARC 
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GO TO 56 
54 J2^Y = 2 
IOPARC = I 
ISWTCH = J2hY 
I F ( X2WAY . EQ. 1 ) GO TO 53 
56 CONTINUE 

IKPOP2 = PCAPI IOPARC) 

IKPOP1 = PCAPt CPTARC ) 

MCAP ( OPTARC ) = NEWC P ( CPT ARC ) 

KCAP ( IOPARC) = NEWCP ( I CPARC ) 

C STRIKE BEST ARC AND FIND FLOW CCST 

CALL OCK( NUMNDE ,NUPARC,FRGM,TO, JFCST,MCAP,LRNC, FLCW, 
IP I ,NL, INFEAS) 

I F ( INFEAS. E0.1IRETURN 

CALL FLCCSTIFLOW.JFCST ,NUMARC, IMNCIJT , I NVL, 1 ) 

HCAP (OPT ARC ) = IKPCP1 
F'CAP ( I CPARC ) = IKPCP2 
RETURN 
END 
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SUBROUTINE OUT PUT ( J ) 

COMMON/D AT 01 /I NVL ,DAY, NUMDAY, NUMARC, OPTARC,L INE, INFEAS 
C0MM0N/DAT02/S,NUMNDE, MCAPU,K£LOI , NR,N 
COMM ON /OAT 03/NR Ct I MNCST < IMNCUT , X2WAY , MNCO ST, MNTC ST 
COMMON/DAT 0 A/REPDA Y( 150) ,ARNAM (1 50 ) , FROM ( 150 ), TO ( 150 I 
COMMON/ CATC 5/RMAXCT( 150) ,RPMCST( 150) ,PI ( 150) 
C0MMCN/DAT06/LBND( 150 > , NEWCP ( 1 50 ) , MC AP ( 1 50 ) , IS( 150) 
COMMON /O AT 07/MRTM( 150) ,NN(15C) ,MCTY(150> , JFCST(150) 
COMMON/ DAT 0 8/ FLOW ( 150) ,NL( 150) 

COMMON/D AT OS /VAL UE ( 150) , PRONAM(25) ,J FLOW (150) 

COMMON/ DAT 1 0/REQFL0,REDU,J 2 WY, I SWTCH ,CRC 
C0MM0N/DAT11/I0PARC,T0TDST(150 ) t TOTCST( 150) , IT(150) 
COMMON/DAT 12/ 1 V( 150) , DAMAGE (150) , VP (150 ) 
C0MMCN/DAT13/IRQFL0,VULCAP( 150) , CRT ( 150) 

INTEGER TOTCST,TOTCST ,DAY,ARNAM,X2WAY, FROM, TO 
INTEGER FLOW, PI , OP TARC , REPDAY , S , RMAXCT,CRC 

IF( JoGTol ) GO TO 2C 

WRITE (6, 506) ( PRONAM ( I > ,1=1 ,20) 

I F ( REDUo LT© lo ) GO TO 10 
WRIT E ( 6 , 507 ) NUMDAY , N 
LINE = 8 
RETURN 

10 WRITE ( 6, 510 )N, NUMDAY, FL0W(1 ) , I ROFLO, I ROFLO , N 
WRITE(6, 50 A) 

LINE = 0 
RETURN 

20 I F(L INE 0 NE 0 A 8 ) GO TO 30 
WRITE ( 6 , 50 A ) 

LINE = 0 

30 I F ( L INE 0 NE 0 OGO TO *0 

AO GO T0(10 0,100, 150,160, 210, 50, 1 80 )»J 

50 I ARC = 0 

I F( DAYo EQoO )WRITE( 6,502) 

I F( DAYo EQ 0 0 )G0 TO 70 
WRITE(6,519)DAY 
60 LINE = 0 
70 WRITE(6,508) 

80 I ARC - IARC + 1 

IF(IARCoGToNUMARC) GO TO 90 
IF(FLOW( IARO 0 EQ 0 O 0 AND 0 DAY 0 NE 0 O) GO TO 80 
LINE = LINE + 1 
II = FROM ( IARC) 

I J = TO ( IARC) 

WRITEt 6, 509) ARNAM( I ARC) ,NN< I I ) , NN( I J ) ,LBND( I ARC) , 

1 MCAP ( I ARC) , JFCST(IARC) ,RMAXCT< I ARC ) , MRTM ( I ARC ) , 

2FL0W ( I ARC ) , VP( IARC ) 

I F ( L INEo NEo A9 ) GO TO 80 
WRI T E( 6 , 50 A ) 

GO TO 60 
90 CONTINUE 
LINE = 0 
RETURN 

100 CONTINUE 

I F ( L I NE + 60 GEo A 8 ) GO TO 200 
IF ( L INEo NEo 0 ) GO TO 120 
1)0 CONTINUE 

WRITE(6,5C3)DAY 
WR I T E ( 6 , 500) 

LINE = LINE + A 
120 CONTINUE 
K = OPT A RC 
III =FR0M ( K ) 

J JJ=TO( K) 

WRITE< 6, 501 )NR,NN( I II ) ,NN( JJJ) , MCAPU , MC AP ( K) , 

1JFL0W( K) ,Kf LOI, REPDAY ( K ) , CRC , FLCW( 1) , IMNCUT 
2 , 1 S ( K ) 

130 LINE = LINE + 1 
RETURN 

150 WRITE(6»513)DAY 
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LINE =L INE + 8 
RETURN 

160 WRITF! 6, 512 )DAY 
LINE = LINE + 8 
RETURN 

180 WRITE! 6 , 51 5 ) IMNCST , FL 0W ( 1 ) 

I F ( DAY o EOo 0 ) GO TO 190 
IF(FL0W(1 )oGT o RE0FL0)G0 TO 190 
WRITE! 6, 51 6 ) DAY 
190 CONTINUE 

WRITE! 6,504) 

LINE = 0 
RETURN 

200 L1NF =48 
GO TO 20 
210 CONTINUE 

WRITE(6,520)DAY 
WRITE (6,521) DAY 
DO 220 J ARC=S , NUMARC 
KK = FROM! J ARC ) 

KJ = TO(JARC) 

IFITOTDST! JARC)o EQoOo ANDoTOTCST ( JARC ) o EQ o 0 ) 

1 GO TO 220 

I F ( F LOW (JARC)oEOoO) GO TO 220 

WRITE! 6,522 )ARN AM! J ARC) ,NN( KK) , NN ( KJ ) , T OTDST ( J ARC ) , 
1TOTCST ( JARC ) 

220 CONTINUE 

WRITE(6,523)NR,NRC 

RETURN 

500 FORMAT (T30, 'ARC CA PAC I TY • , T48 , • ARC FLOW • , T 61 , • TO BE*, 
IT 71 , 'CUMo * ,T77, • MAXIMUM* ,T87 , • NETWORK* ,T96 , • NUMBER* , / 
2, T10, 'MISSION* ,T28 , 'BEFORE AFTER ', 4X ,* B EFOR E AFTER' 
3, 2X, 'REPAIRED* ,2 X, 'REPAIR* , IX , • NETWORK* ,4X, ' FLOW* ,4X, 

4 * OF ARC* ,/ ,T10, * NUMBER ', 2X, 'FROM TO' ,2X, ' STRIKE' ,2X, 
5* STR IKE* ,3X , • STP IKE STR IK E • , 2X , • ON DAY ', 4X, • COST *, 3X , 
6 'FLOW* ,6X, 'COST* , 4X , • STR I KES • ) 

501 FORMAT (10X, 1 3 ,5X , A 2 , 3X , A2, 2X, 2 ( I 5» 3X) , 1 X, 2 ( I 5, 3X) , 

114, 5X, 15, IX, 17, IX, 19, 3X, 14) 

502 FORMAT! *0' ,19X, • INITIAL FLOW PATTERN FOR LOC NETWORK* 
It/) 

503 FORMAT ( *0* ,8X, 'STRIKES TARGETED FOR DAY', 13, 'o') 

504 F0RMAT(1H1,10X) 

506 FORMAT ( 1H , 28 X,20A4 ,/// ) 

507 FORMAT ( 1H ,32X, 23HCAMPAIGN WILL LAST FOR ,13, 

112H DAYS, WITH ,13,' MISSIONS PER DAYo * ) 

508 FORMAT ( 44X , 'LOWER* ,6X, 'UPPER' , 13X , • MAXo • , 5X, 'MAXo • , 1 5X 
1 , ' ARC' , / , 20X, ' ARC' ,4X, • FROM' ,5X,*T0*,6X, 'FLOW' ,7X, 

2' FLOW* , 5X, 'FLOW' ,4X, 'REPAIR* ,3 X, 'REPAIR* , 4X , ' ARC * ,7X, 

3 * VULNERo ' , / ,20X, 'NAME' ,3X, 'NODE* ,4X, 'NODE' ,5X, 'BOUND* , 
46X , ' BOUND* ,4X, 'COST • ,5X , 'COST' , 5X, 'TIME* , 5X, ' FLOW* , 

56X , 'PARAMo ' ) 

509 FORMAT ( 2H , 1 8 X, A4 ,4X , A2 , 6X , A2 , 1 X , I 8 , 5X, I 8 , 1 X, 1 7 , 

1 2 ( 4X , I 5 ) ,1X,I8,6X,F5o3) 

510 FORMAT ( 36X, • MAXIMUM NUMBER OF STRIKES EACH DAY :• 

1, 13, /,36X, 'NUMBER OF DAYS IN THE CAMPAIGN : • ,1 3 , 

2//,36X, 'NOTE: MISSIONS WILL CONTINUE EACH DAY UNTIL', 
3/ ,36X, ' ONE OF THE FOLLOWING OCCURS :',/,36X, 

4 ' ( 1 ) UNINTERDICTED FLOW (',14,' UNITS) IS',/,40X, 

5* REDUCED TO ',14,' UNI TSo ' , / , 36X, • ( 2 ) FLOW AT THE ', 
6'BEGINNING OF E ACH • , / , 40 X, ' DAY IS REDUCED TO ',14, 

7* UNITSo ' , /,36X, • ( 3) ',13,' STRIKES HAVE BEEN ' 

8, ' TARGETEDo • ) 
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512 FORMAT ( 1H , / / / 10 X , 2 7HF LOW TOTALLY STOPPED ON DAY , 
1110,/// ) 

513 FORMAT ( 1H , ///10X, 30HMIN IMUM THRU-PUT CANNOT BE MET, 

1 • ON DAY* ,110,///) 

515 FORMAT (' O' *// ,20X, ' FLOW COST :',I8,/,20X, 

1 * THROUGHPUT :» ,18) 

516 FORM AT ( * 0* , / / , 20 X, ’NO STRIKES REQUIRED ON DAY*,I3,'o*» 

519 FORMAT! ' O' , 19X, • ALL ARC REPAIR SCHEDULED FOR DAY*, 13, 
1* HAS BEEN COMPLETEDo ' , / , 20 X , • RESULT I NG FLOW PATTERN* 
2, * IS AS FOLLOWS: • , /) 

520 FORMAT! • O' , /,20X, • RECAP FOR DAY ', I 3 ,':»,/, 20 X , 

1 • #*******♦****:*••♦:** • j 

521 FORMAT (• 0* ,45X, • TOTAL NUMBER TOTAL NUMBER ',/, 2 1 X , 

1* ARC* , 21X , * ARC STRIKES ARC STR I KE S * , / , 2 2X , • N AME ’ , 

24X, ' FROM* , 5X, » TO* , 5X, ' FOR DAY* , I 3, 6X , • FOR CAMPAIGN') 

522 FORM AT ( 2H , 20X, A4 , 5X , A2 , 6X , A2 , 9X, 1 4 , 1 2X , I 4) 

523 FORMAT! *0* , /,22X, * NUMBER STRIKES THIS DAY :*, 

119, /,22X, 'NUMBER STRIKES THIS CAMPAIGN : * ,19) 

END 
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APPENDIX D 



COMPUTER OUTPUT FOR SAMPLE PROBLEM 
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